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Abstract 

Here, cosmology is obtained from the variable gravitational constant 
G oc 0~ 2 with 4>{x) being a scalar and its fluctuations around the 
ground state. The gravitational action contains Einstein-Hilbert like 
term with variable G, kinetic energy and self-interaction potential for 
4>{x). Two phase transitions take place in this model. The first one 
takes place at the GUT (grand unified theory) scale ~ 2.45 x 10 14 GeV, 
when the early universe exits the inflationay phase and the second one 
at the electro-weak scale. Spontaneous symmetry breaking takes place 
around this scale As a consequence, variable G acquires constant value 
Gn (the Newtonian gravitational constant). The standard model of cos- 
mology is obtained in the post-second phase transition era. Interest- 
ingly, the dark matter and quintessence dark energy are created from 
the gravitational sector as a combined effect of the linear term of scalar 

curvature and (f)(x) without using non-linear terms of curvature or any 

l 
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other matter fields. Radiation-dominance is obtained in both the pre- 
second phase transition era and the post-second phase transition era. 
In the pre-second phase transition era after inflation, radiation domi- 
nated phase of the universe decelerates in the non-standard manner as 
~ t 1 / 3 , but expands in the standard way as ~ t 1 / 2 in the post-second 
phase transition era. It is interesting to find that fluctuations in <j)(x) 
cause the phantom dark energy after the second phase transition, which 
grows later in the expanding universe. Dominance of dark matter over 
radiation is obtained at red-shift 1091.89. Matter-dominated phase 
decelerates as ~ t 2 ^. Phantom dark energy, created during radiation- 
dominance of the post-second phase transition era, overpowers cosmic 
dynamics at red-shift 0.46. The late universe accelerates driven by the 
phantom dark energy, when time t > 0.13£o (where to is the present age 
of the universe). It is found that phantom dark energy evaporates ul- 
timately contributing to the vacumm energy with density almost twice 
of the present dark energy density. It is also found that the phantom 
era is free from big-rip singularity. 

Key words: Variable G, quintessence dark energy, early inflation, phase 

transitions and phantom dark energy. 

PACS nos. :98.80.-k, 04.50+h, 98.80. Cq, 98.90.+S 

1. Introduction 

Observation of the late cosmic acceleration [1, 2], made during last 
few years, challenged the standard model of cosmology predicting decel- 
eration in the present universe. This revolutionary result (late cosmic 
acceleration) prompted cosmologists to think dominance of some exotic 
fluid in the very late universe. Theoretically, the basic characteristic 
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of such a fluid is obtained from Freidmann equations giving the cos- 
mic dynamics. These equations suggest negative pressure p < — (l/3)p 
with p being the energy density) for the dominating fluid to get the late 
cosmic acceleration. This requirement implies violation of the strong 
energy condition (SEC). Being unknown and having gravitational ef- 
fect, this exotic fluid is popularly known as dark energy (DE). Later on, 
the case p < —p was found better fit to observations. In this case, the 
weak energy condition(WEC) is also violated. This type of dark energy 
fluid was dubbed phantom [3]. Apart from these models, many other 
field-theoretical and hydro-dynamical models of DE were proposed in 
the recent past to explain challenging observation of late acceleration. 
A review of these investigations is available in [4]. In the race to search 
models to explain the phenomenon of late cosmic acceleration as well 
as to satisfy solar system constraints, curvature was proposed as gravi- 
tational alternative of DE. This work is reviewed in [5]. These models, 
specially model with non-linear terms R n and R~ m [6] received criticism 
in [7] on the ground that these models did not produce matter in the 
late universe needed for formation of large scale structure in the uni- 
verse. In another review, Nojiri and Odintsov discussed dark matter 
in the late universe responding to criticisms against their work [8]. In 
refs. [5, 7, 8], non-linear curvature terms are treated as DE term a 
priori and its consequences are discussed. But, instead of assuming 
curvature as DE term it is more appropriate to know whether curva- 
ture can produce DE density terms in the Friedmann equation without 
assuming non-linear curvature terms as DE term a priori. This aspect 
is addressed to in [9, 10, 11, 12, 13, 14] using an approach different from 
the approach of refs. [5, 7, 8]. In this approach, matter is created from 
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gravity itself, which resolves the problem raised in [7] against modified 
gravity . 

In spite of so many efforts to investigate an appropriate model of 
DE, it is not unreasonable to think for some other alternative, which is 
planned to probe in what follows. In this context, it is useful to consider 
our past experiences. We know that Kaluza-Klein (K-K) modes are not 
excited below a TeV scale. According to Ref. [15], string/M-theory 
excitations are also expected at a TeV scale. Moreover, coupling of 
K-K modes to matter or gravitational excitations is set by the weak 
scale. Electromagnetism and weak interaction are unified at the scale 
246 GeV, which is observed. Although at CERN's LHC it is possible 
to reach the energy scale 10 4 GeV, we have a renormalizable theory at 
246 GeV in the observed sector (GUT scale is not observed) of mass 
scale. 

So far, gravity is observed much below this scale. One possible rea- 
son for non-observation of gravity at high energy scale is the effect 
of different non-gravitaional forces at high energy, which are almost 
negligible at low energy. 

These arguments, which are based on experimental and theoretical 
results, prompt us to think deviation in the gravitational constant G = 
Gn ■ One possible deviation is the space-time dependence of G. On the 
basis of above arguments, this type of deviation is proposed here above 
the electroweak scale 246 GeV. Another motivation for variable G is 
given by the possibility of physical constants to have different values in 
the remote past of the cosmic history. An evidence for this possibility 
is given by the apparent observation that fine structure constant had 
different value in the distant past [16]. 
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This type of proposal is not new in the history of gravity. Varying G 
had been considered by many authors in the past [17, 18, 19, 20, 21, 22]. 
First time, it was proposed in Ref.[17] by Jordon in his scalar-tensor 
theory taking G proportional to an arbitrary power of 0" 1 with <p(x, t) 
being a scalar. This theory was developed later by Brans and Dicke 
in Ref. [18] using G oc <p~ l . In these theories, the kinetic energy term 
for <f)(x,t) is non-canonical having the form u;<9 M 0<9 M 0/0 with oo being 
a coupling constant. Later on, in Refs.[19, 20, 21, 22], these theories 
were generalized considering oj dependent on 0. 

In the old gravitational theories with variable G (OGTVG) proposed 
in the past [17, 18, 19, 20, 21, 22], G is not equal to Gn even in the 
present universe. But experiments support G = Gn in the present 
universe as results of Einstein's theory of gravity, a theory with G = 
Gn, support observations of advance of perihelia and bending of light 
rays around the Sun satisfactorily , whereas the Brans-Dicke theory 
which is the most popular theory having variable G in the present 
universe too has problem in explaining advance of perihelia and bending 
of light rays. Astrophysical observations, made at the dawn of this 
century [1, 2] concluding late cosmic acceleration as well as content 
of the present universe in the form of 73% dark energy and 23% dark 
energy, are based on G = Gn in the late universe. 

Moreover, on taking G ~ the kinetic energy term for is non- 
canonical. So, a theory is needed to explain variable G in the past and 
G = Gn in the present universe. 

It is possible to have a theory with variable G at high energy scales 
but G = Gn in the late universe and canonical kinetic energy term for 
(f) , if G ~ is replaced by G ~ 0~ 2 with potential V(4>) = A(0)|0 2 — 
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/3 _1 Mj>| 2 /4, where dynamical coupling A oc 0~ 2 is dimensionless, Mp 
is the Planck mass and f3 is a dimensionless coupling constant. This is 
possible without causing any harm to the theory. 

On the basis of above experimental probes getting a renormalizable 
theory upto the electro-weak scale 246 Gev in the particle sector and 
probe of gravity upto 0.1mm., it is proposed here that phase transition 
takes place at = p~ 1/2 \M P \ giving G = f3<p- 2 = Mp 2 = G N at the 
electro-weak scale. But prior to the phase transition, G = j3(j)~ 2 ie. 
above the energy scale 246 Gev, G is variable. 

The present model is different from OGTVG proposed in the past 
[17, 18, 19, 20, 21, 22]in the following four different ways. 

(i) G oc 0~ 2 in the present model, whereas G oc 0" 1 in OGTVG except 
Jordon's theory with an arbitrary power of _1 . 

(ii) In the present model, has canonical kinetic energy term and quar- 
tic self-interaction potential with a dynamical coupling. This set-up is 
in contrast to OGTVG having non-canonical kinetic energy term and 
no potential for 0. 

(hi) In the present model, G is a variable above the electro-weak scale 
and this is equal to Gn in the present universe. This set-up is in 
contrast to OGTVG , where G is always a variable, 
(iv) Mass dimension for is onein the present model, whereas mass 
dimension of is two in OGTVG (in natural units given below). 
Set up of the model and consequences. Here, the gravitational 
action contains Einstein-Hilbert like term with G = (3(j)~ 2 (with (3 being 
a dimensionless coupling constant), kinetic energy for (required for 
to be a dynamical variable) and self- interaction potential V(<p). is 
extremely small at the Planck scale and increases with the expansion 
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of the universe. It is obtained that, in the flat, homogeneous and 
isotropic model of the universe, trace of gravitational field equations 
derived from the considered action yields the Friedmann equation with 
dark matter density and energy density of an unknown fluid, which 
mimics quintessence dark energy and derives early inflation. These 
energy density terms emerge from the gravitational sector ( mentioned 
above) without taking any other source except and scalar curvature 
R. 

Using observational consequences, it is found that universe exits from 
the early inflation at = C = 3.45 x 10~ 13 M P expanding around 10 27 
times. 

Two phase transitions take place in this model. The The first one 
takes place at a = a c (when <p = <f> c ) and enormous amount energy is re- 
leased as latent heat. The emitted radiation (latent heat) heats up the 
universe upto the temperature T c = 2.45 x 10 14 GeV. Due to creation 
of this radiation , the derived Friedmann equation gets modified. It 
is found that the emitted radiation dominates the effective Freidmann 
equation after inflation. As a consequence, the universe decelerates 
with a(t) ~ t 1 / 3 , which is different from the standard radiation model 
due to the effect of variation in G. 

The second phase transition takes place at = (3~ l ^ 2 Mp and energy 
is released again. It is found that this phase transition (the second 
one) is a bit different from the first one. At the epoch of the second 
phase transition, major portion of the energy is released as radiation 
heating the universe slightly by 5.22 x lO^fr above the temperature 
246GeV, which is the temperature of the radiation being red-shifted 
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from the epoch of the first phase transition. An extremely small por- 
tion of the released energy causes small fluctuations in the stable state 
4> = fi~ l l 2 M P . These fluctuations die out soon causing creation of 
the phantom dark energy. Radiation still dominates dynamics of the 
universe but with a change during the post-second phase transition 
era. In this era, radiation phase expands as a(t) ~ t 1 ' 2 consistent with 
the standard radiation model. It is because spontaneous symmetry 
breaking takes place at the second phase transition and G becomes a 
constant equal to Gn = Mp 2 . Thus the standard cosmology is recov- 
ered at this epoch (the epoch of the second phase transition) with the 
background radiation having the initial temperature ~ 246GeV. 

This scenario is followed by matter-dominated universe expanding 
with a(t) ~ t 2 / 3 . The created phantom dark energy dominates the 
universe at red-shift z = 0.46. As a consequence, universe speeds up 
due to a jerk caused by phantom energy. 

Natural units H = c = kp = 1 (where H = h/2n with h being 
the Planck's constant, c being the speed of light and ks being the 
Boltzmann's constant) are used with GeV as a fundamental unit. In 
these units, it is obtained that lGeV = 4.96 x 10 13 °ir = 9.07 x lfr 24 gm 
and lGeV" 1 = 3.879 x 10- 15 cm = 1.293 x 10- 25 sec. Here, Newtonian 
gravitational constant is given as G N = Mp 2 with M P = 2.4 x 10 18 GeV 
being the Planck mass. 



2. Early inflation and pre-second phase transition era 



Investigations begin from the action 
S — d^xsf^g 



16tt 2 y ^ vV ! 
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where R is the Ricci scalar curvature , 

A(0) = K<j>- 2 (lb) 

with constant K having mass dimension 2 and <9 M = d/dx^. 

Using the condition SS/Sg^ = 0, action (la) yields gravitational 
equations 

] 



2 1 



16tt' 



A, 



with 



□ 



-5 



^<rd v 



(2b) 



Taking trace of (2a) and setting (3 = 8n/3 for convenience of calcu- 
lations, it is obtained that 



~R - [<f>U(f> + d^(pd u 4>] + \d^<pd u <p + ^ 



6 ~ LT T ■ r "^ J ' 2" ' 2 7 8tt ^ °' (3) 

Moreover, using the condition 5S/5<p = 0, action (la) yields equation 
for <fi as 



□0 + ^i? + A0 



+ 2A0" 1 



= 0. 



(4) 



Eliminating D<f> from (3) and (4) as well as re-arranging terms, it is 
obtained that 
1 



-R-cfi-'dyd^ = 2A0- 



~-4> 2 




5 


2 - 




2- 






~ 2 







(5) 



Experimental evidences support spatially homogeneous flat model of 
the universe [24]. So, the line-element, giving geometry of the universe, 
is taken as 

dS 2 = dt 2 - a 2 (t) [dx 2 + dy 2 + dz 2 } (6) 
with a(t) as the scale factor. 
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In this space-time, the scalar curvature is obtained as 



R = 6 



a /a\'- 
a \aJ 



(7) 



Due to spatial homogeneity of the space-time (6), <f>(x,t) = <f>(t) and 
A(0) = X(t). As mentioned above, in the proposed model 0(t) increases 
with the expansion of the universe. So, it is reasonable to take solution 
of (4) as 

«*> = {iMd* (8a) 

with n being a positive real number and a* being the scale factor when 
= \J^\Mp\ an d phase transition takes place, which is discussed 
below. 

Moreover, A(i ) = A is proposed to be 

A(t) = A p) 2n , (86) 

where ao = a(to) with t being the present time. 1(b), (8a) and (8b) 
yield 

^ = ^AoM|(^) 2 ", (8c) 
Connecting (5), (7), (8a) and (8b), we obtain 



a fa\ 



+ - 

a \a/ 



—n 



2K 



8tt) + 27rla 



2 la 



(9) 



First integration of (9) yields the Friedmann equation (FE) 
'd\ 2 C _/a*\ 4 «r 1 / 3 \ 2 



„<-. +2A "(7 



(a + 2 -An) \8ir 



2n(a + 2 - 2n) V a 



a* 



2rj 



2(a + 2) V a 



a* 



An 



(10a) 



where, 

a = 2-n 2 . (106) 

In FE (10a), the term C/a^ 2+a ^ corresponds to matter density for 
a — 1. So, to get the matter , which is required for formation of large 
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structure (planned to be discussed in other paper) in a viable cosmology 
a — 1 is taken. As a result, 

n = 1 (10c) 

from (10b). 

Further, using n from (10c), (8a), (8b) and (10a) are re- written as 



7r \a 



A(t) = A (^) 2 (life) 

and 

Thus, it is interesting to see that the matter term C/a 3 in (12)emerges 
spontaneously from the gravitational sector, so it is recognized as dark 
matter (DM) being non-baryonic. The spontaneous emergence of dark 
matter is also obtained in references [9, 10, 11, 12, 13] from higher-order 
terms of scalar curvature. The present DM density of the universe is 
observed as p m = 0.23p° r [25] with 

pi = M = 2.85 x 10~ 50 GeV 4 (13a) 
8ttG 

and 

H = — = 2.035 x 10~ 43 GeV (136) 
to 

with t = 13.7Gyr = 6.6 x 10 42 GeV _1 and G = Mp2,M P = 2.4 x 
10 18 GeV in the present time. 

Dark matter is characterized as pressureless fluid, so energy conser- 
vation equation 

p + 3-p = (14) 

CI 

yields 

0.69# 2 Mj,/ao\ 3 
Pm = Svr UJ (15a) 
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using observational values given above. The DM density obtained 
above is 

Pm= w=^U (156) 

using (11a). At a = a*, we obtain 

9M P C _ 0.69# 2 M 2 /a 



*)* (15c) 
a*/ 



(87r) 2 a3 8tt 
from (15a) and (15b). (15c) yields the integration constant C as 

C = 184 f°°° 3 . (15d) 

Like the DM term, other terms in (12) given as 

also emerge from the gravitational sector. Being caused from the grav- 
itational sector, pc is recognized as density for the dark energy fluid 

According to (16), pc = at a = a** given by 



1 (a**\ //l /3 



y/2 \ a, J V Vtt V 2/ 27TV/2 

Using the pressure density p = up with being the equation of 
state (EOS) parameter in (14), EOS parameter for the dark energy 
with density pc is obtained as 

2 . • (16c) 



2 4a 2nV2 



3 3d* / a* \ 



(°a) _(^\/^) {v^(«*) W2} 



(16c) implies ujq < —1/3 for a < a*. It means that this dark energy 
mimics quintessence. 
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In the very early universe, 
C _ 1.847r# 2 ^a y 

Using the inequality (17), (12) is obtained as 

/ (X \ ^ / CL^ \ ^ / fl c \ ^ / 1 / 3 \ ^ fl / CL \^ 3 1 ^ / ^* \ ^ 1 ^" 

W ~ 2K W/ V~a~J Lw V 2/ "l^Va^J "2^7^/ 1^/ J J' 

(18) 

where a c is the scale factor at which early inflation ends. 
(18) integrates to 

sechO - tanhQ = De 2 ^, (19a) 

where 

i fi ccche _ 1 ( a Y 3 } 2 ( a *Y am 

and 

Va c / 87r Va*/ Va c / 

(19a) and (19b)yield 

/ay _ ijVM» = F/apy UV^l^) {19d) 

\a c J 2ix) \a c J l\a c J 2ttv / 2- 1 J 

with fT c given by (19c). 

For dp < a < a c « a*, (19d) is approximated to 

a ~ a P e^ (i - ip) , (20) 

showing the early inflation. 

As a << a*, the slow-roll parameter e is obtained as 



' = -£4 (21) 



where if = ^ is given by (18). 
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In Refs. [27, 28], formula for maximum number of e-fold during 
inflationary phase was obtained as 



63.3 + 0.25/ne ~ 63.13. 



Connecting (21) and (22a), it is obtained that 

a c 



~ e 



63.3 



3.1 X 10 



27 



dp 



(22a 



(226) 



Exit from inflation and radiation-dominance 



As mentioned above, inflation ends at time t — t c ,(f> — (f) c and the 
first phase transition takes place. At this epoch, an enormous amount 
of energy is released as latent heat with density 

.2 /n n \2 3 2 

87 



v{4> P ) - v(<f> c ) 



\ M P 



)l( 



/ao y (^p\ _ 
\apJ \ a* / 

2 3 2 n 



8tt 



4 ap 



ap\' 



(23) 



According to WMAP results [25], density of the released energy is 
obtained as 

V{4> P ) - (0 C ) = — x 10- 18 M P = 3.96 x 10 54 GeV 4 (24) 
8n 

at = 4> c . 

The emitted radiation (latent heat) thermalizes the universe after 
inflation upto the temperature T c = 2.45 x 10 14 GeV, which is the GUT 
(grand unified theory) scale. T c is given by 



ILr 4 

15 c 



= 1.19 x 10 57 GeV 



(25) 
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using (23)and (24). 

Thus energy with radiation density, given by (25), is emitted at the 
end of inflation and 

At the end of inflation, the density of the left - over 0— energy is 
obtained as 

V(<p c )= Prem = ^A Mp— y 
4 a c J 

15 c \aJ 
= 0.412GeV 4 (26) 

using (24) and (25). 

This energy is released after the second phase transition at a = a* 
being discussed in the following section. 

The emitted radiation with temperature T c = 2.45 x 10 14 GeV (eval- 
uated below) is given by density 

A-£?(f) 4 - <w.) 

Obviously, contribution of this radiation affects the cosmic dynamics 
leading to modification of the FE (12) as 

( 1 / Qj \ ^ 3 *\ 2 x x 2 *\ 2 - 

Using (11a), (lib), (25) and (27a), it is obtained that at a = a c , 

^>-© 2 e) 4 [(^) 2 -{^© 2 -^} 2 © 2 } 2 ]- 

(27c) 

Using inequalities (17) and (27c), (27b) reduces to 

/ay ^_8^ 4^ md) 
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incorporating (lib) and (27a). This equation is integrated to 



T /87T 3 T 2 a2 a T 1/3 

a(t) = a c [l + -^^fH - t c )\ (29) 

The result (29) justifies exit of the early universe from inflationary 
phase at t = t c . It also shows deceleration for t > t c . 

Here, it is important to remark that the expansion (29) is different 
from the standard radiation model due to the effect of variable G = 

3. Post-second phase transition era 

It is discussed above that the first phase transition takes place at 
= C . But this state is not stable. So, rolls down the potential hill 
V((f)) to a stable state being obtained from the condition dV {(p) / d(p = 
. For the potential V(4>) in action (1), this condition looks like 



3 



8tt 

(30) yields stable fixed points 



Mp) <f) 2 -^-M' P =0 (30) 



8vr 



= ±J£-Mp. (31) 

V 07T 

Spontaneous symmetry breaking takes place, when <p reaches the 
stable states given by (31). As a consequence, G acquires a constant 
value Mp 2 = G N (the Newtonian gravitational constant). Moreover, 
the second phase transition takes place and energy with density V((j) c ) 
is released at = yJ^\M P \ . In what follows, it is discussed that the 
major portion of this energy is released as radiation, but an extremely 
small portion causes fluctuations in the state = ^J^\M P \ dying out 
in a very short interval of time. 

So, is taken as 

= J ^-M P + 50, (32) 

V o7T 



VARYING GRAVITATIONAL CONSTANT, 



17 



where 5(p is a small fluctuation in the state = 
Using (32) in (3) and (4), it is obtained that 



\M P 



{M\ + J^M P 8<t>) 



R+\\l -M P n5(f) = 

2 V vr 



and 



aS(j) + -R5<j> + 2K5(f) = 



(33) 



(34) 



neglecting higher orders of 5(f). 

In the space-time (6), (34) is re- written as 
a 



50 + 3- 



+ 



rd_fa 
J \dt\a 



-J) 



+ 2K 



0. 



(35) 



(35) is an equation of damped harmonic oscillator having decreasing 
amplitude with increasing scale factor a(t), so a trial solution can be 
taken as 

S(P = - (^Jjcos[w(t - U)} (36) 

where r is a positive real number, w is the frequency and t* is the phase 
transition time. Here a e « a is the scale factor, when fluctuation 
vanishes. 

(36) satisfies (35) if 

3 IK 



r = 



2' 



zu = 



with a being a constant and 

(!) 



5 /aN2 
a/ 2Vct/ 



2{w 2 -2K). 



Connecting (7), (33) and (36), we obtain 

2 2w 2 12k 



a /a\- 
-+(-) 




(37a, 6) 



(37c) 



(38) 



after some manipulations. 
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(37c) and (38) yield 
4 



a 



a 



9 n T ^ Aw 1 

-(ro 2 - 2X) 

3 V ; 3M P 




4 " 
3 . 
AK 
3a 



2 - 2AT) - 
1 - 2a) - 



2tt 
2~^ 




T 

2tt 




a 

a 



a*\ 3 / 2 y 



a*\ 3 / 2 



using (8c), (10c), (36) and (37a,b). 

The kinetic energy density of the oscillating 



is obtained as 



1 

PS4> = 2 [ 



— ^— ) + zutan[zu(t — t*)] 



(40) 



using (36) and \cos[w{t — t*)]\ < 1. 

(40) shows that kinetic energy of 5(f> decreases with growing a(t). It 
means that during oscillation, 5(f) emits energy, which is being pumped 
to the universe. The energy density in the universe caused by fluctu- 

J^|Mp| is given by the energy 



ations around the stable state <j) 
density 



Pphde 



KMp 
2an 



(1 -2a) - 



Mr 




3 I V a 



}] (41) 



As phantom dark energy is seeded by decay of 5(f), so p p hde = at 
a = a*, p p hde = at a = a*. Using this condition in (41) 



(l-2a) = 



A 2n 




M P V 3 I 

Now, (8c), (10c), (41) and (42) yield 



a*\ 3/2 



(42) 
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When Pphde, given by (43), satisfies the conservation equation (14), 
the equation of state parameter w p hde is obtained as 

3/2 



w phde = -1 3/2 • (44a) 

2[l-0) 

The result (44a) implies that the energy emitted due to decay of <50 
mimics phantom dark energy as w p hde < — 1- 

Further, according to observations, dark energy density of the present 
universe p phde = 0.73p° r is extremely low and phantom energy density 
increases with the expansion of the universe. So, phantom energy den- 
sity is expected much lower than p phde for a < a . It shows that an 
extremely small portion of released energy with density 

e = 10- x V((f) c ) = 0.412 x lO^GeV 4 (446) 

is converted into phantom dark energy and the major portion « 0.412GeV 4 
is emitted as radiation at the second phase transition like latent heat. 
Here x is a positive real number. 

As discussed above, the emitted radiation after the first phase tran- 
sition has temperature 246GeV after red-shifting to the scale factor 
a* > a c . Using these value in the formula for temperature in the ex- 
panding universe, we obtain 

- = £- = 10 12 . (45a) 
a c 246 v ; 

Connecting (8a) and (45a), it is obtained that 

— x l(T 12 Mp = 3.45 x 10~ 13 M P . (456) 
8tt 

Moreover, at the second phase transition, radiation is emitted hav- 
ing density given by (26). Now, at this epoch, radiation is available 



20 S.K.SRIVASTAVA 

from two sources (i) red-shifted radiation emitted at first phase tran- 
sition having temperature 246GeV and (ii) radiation being emitted at 
the second phase transition. So, total radiation density at a = a* is 
obtained as 

„* _ 71-2 rp4 

Prd " 15 i* 
7T 2 

= — x (246) 4 + 0.412. (46a) 

(46a) yields 

T* = 246 [1 + 4.275 x 10~ n ]GeV = 246GeV + 5.22 x 10 5 °K. (466) 

This result shows that , after the second phase transition, the universe 
is reheated slightly by temperature 5.22 x 10 5 °K. 
So, in this era, radiation density is obtained as 



As the released energy causes fluctuation in the state 4> = \l ^-\M P \, 



8(f) has energy with density 

p 54> = 0.412 x lO^GeV 4 (47) 
at a = a*. This result is obtained using (44b). 

Radiation-dominance during post-second phase transition era 

In this era too, radiation-dominated universe is obtained ,but with a 
changed scenario. In this case, Friedmann equation (27d) gets modified 
again and the variable G is replaced by the constant G = G N = M p 2 
leading to the standard form of FE during radiation-dominance. The 
modified (27d) is obtained as 



8tt 3 _ T4 fa*\ 4 



AhMl * V a 



= .i.M. : l()- 27/ '" 
a 
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using (25), (45a) and (46c). 
(48) is integrated to 

a = a*[l + 1.19x 10~ 13 (t-^)] 1/2 (49) 

It shows that the standard model of cosmology is recovered after 
the second phase transition, so emitted radiation with temperature 
T* = 246GeV + 5.22 x lO 50 ^ is identified as the cosmic microwave 
background (CMB) radiation having temperature T = 2.73°K in the 
present universe. 

Using T and (46b), it is obtained that 



4.47 x 10 15 (50a) 
a* Tq 



(45a) and (50a) yield 

^ = ^ = 4.47 x 10 27 (506) 

as well as (26) and (50b) imply 

A = 7.44 x 10~ 127 . (50c) 

After knowing A from (50c) as well as connecting (19c), (20), (22b), 
(45a) and (50a), inflationary period (in the pre-second phase transition 
era) is evaluated as 

(t c - t P ) = 4.41 x 10 7 GeV" 1 = 5.7 x 10~ 18 sec. (51) 

It is discussed above that extremely portion of the energy having 
density g of 5<p decays to phantom dark energy, so 

-Q = Pphde- (52a) 

As phantom is seeded by decay of 5(f), so at a = a*, p p hde — . Using 
this condition as well as (47), (52a) is integrated to 

0.412 x 10-* - p H = Pphde . (526) 
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Using (40), (42), (48) and the condition at a = a*, a is obtained that 

„ N2 2 0.412 x 10~ x 
(l_ 2a ) 2 = -x — .. (53) 

According to (40), = at a = a e , so we obtain 
using (43). 

Connecting (42), (53) and (54), a quadratic equation for (1 — 2a) is 
obtained as 

(1 _ 2a ) 2 _ (1 _ 2a ) + !^£(^) 2 = 0. (55a) 

(55a) yields 

(1 _ 2a) ^ ( ^ = , 53xl0 -, 2 m 

using (50a) and (50c). 

Further (53) and (55b) imply 

10~ x = 2.95 x 10" 52 (55c) 

The present dark energy density is 

P° phde = 0-73p° r = 0.73 x 3 -^l = 2.08 x lO^GeV 4 . (56) 

So, from (43) and (55), it is obtained that 

(43), (53). (54), (55c) and (57) imply 

3/2 



Pphde | J _ ^o*^ 3// | 

2.08 X 10~ 50 



0.412 x 10"* 



= 171.137. (58) 
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Using (50a) for a*/ao in (58), it is obtained that 



\a e ) 



1 - 



a*\ 3 / 2 



(59) shows that 



Thus 



171.137 
.994. 



a e = 1.004a* 



\a J 



p phde = 2.08 x 10 



-50 



MO 



o-(sn 

= 2.026 xl0- 50 {l-(^) 3/2 }. 

using (56), (58) and (59b). 
This result yields 



Pphde Pphde 



for a < a and 



Pphde — Pphde 

for a > a using (50a), (56), (59a) and (59c). 
(44a) and (50a) imply 



w 



phde 



_! fe) 



3/2 



1 _ 2. 

a 



3/2-1 



-24 



23 



(59a) 
(596) 



(59c) 



(59d) 



l + 0.9747{l- (^) 3/2 }j (59e) 



(60) 



= -1 - 3.335 x 10" 

This equation shows that EOS parameter Wp hde for phantom energy 
very much closed to —1 in the present universe. 

Using (44a), (59d) and (59e), w p y e and p p hde are tabulated below. 

Table no.l 
w and Pphde are tabulated for increasing a(t) 
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a^/a 


w 


Pphde 


4.45 x 10 ia 


-83.83 


P^e/171-137 


2.8 x 10 15 


-2 


A^fe/2.05 


10 ib 


-1.06 


Pphde/ 1-15 


1091.98 


-1 - 0.0605 x 10~ 18 


^[1.9747- 0.118 x 10-" 




1.46 


-1 - 2.95 x 10" 24 


p^Jl.9747- 5.9 x 10— j 




1 


-1 - 1.673 x 10~' M 







-1 


l-9747p^ e 



Table no.l shows that w p hde = —1 for a — > oo. This is the EOS 
parameter for the vacuum. It means that phantom dark energy will 
evaporate and contribute to vacuum energy with density 1.9747p° /wfe 



as a — > oo. Moreover, it is interesting to see that p p hde grows very 
fast during radiation-dominance compared to the era of dark matter- 
dominance and phantom-dominance. 

Matter-dominance It is obtained above that, even after the second 
phase transition, radiation dominates the universe. As dark matter 
density p m given by (15b) and (15d)as well as radiation density (46b) 
decreases with growing a(t), a stage is reached at a = ad > a*, such 
that 

ze!w— V = —t 4 (—Y(—Y 

15 * Va^/ 15 * Va / Va^/ 
= 0XA{%)' = 0.66 xl0-°(*)' (61«) 

using (13a). 

Thus, p'd < Pm for a> a d . 

Using (50a) and T* = 246GeV + 5.22 x 10 5 K in (61a), it is obtained 
that 

— = 1092.98. (616) 
®d 

It shows the beginning of induced dark matter-dominance at red-shift 

Zd = ^ - 1 = 1091.98. (61c) 
a d 
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This value of red-shift is very much closed to the red-shift z = 1087^2 
for decoupling of matter from radiation according to WMAP observa- 
tions. 

So, for a > a d , (27b) reduces to 

My = o.W 

Va/ cr 

using (16a) and (16b). 
(62a) yields the solution 

a = a d [l + 0.587# (1 + z d f' 2 {t - t d )f 3 (626) 

showing deceleration. 

It is interesting to see that, universe expands like standard matter- 
dominated model. 

Late acceleration It is discussed above that density of the phan- 
tom dark energy caused by fluctuation around the stable state = 
j^|Mp| grows and matter density decreases with the expansion of 
the universe. So, it is obvious for p p hde to overpower p m on sufficient 
growth of a(t). If this transition takes place at a = a tr , 

3 „r ~./a A 3 / 2 

' a tr / V a t 

using (13a), (16a), (16b) and (59c). (63) yields 



0.6555xl0- 50 (^) 3 = 2.0273xl0- 50 {l-3.33xl0- 24 (^) 3/2 }. (63) 



— ~ 1.46. (64) 

It means that phantom dark energy dominates for red-shift z < 0.46 
and the effective Friedmann equation is obtained as 

, 2 8tt 



fay _ 87T 
\a) ~3M* Pphde 

3.015 x 10 86 {l - (^y ,Z ) (65a) 



using (39), (42), (58) and (59c). 
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On solving (65a), it is obtained that 

a(t) = a*cosh 4/s [27.09 + 2.315 x l(T 43 (t - t tr )\ 

= 3.26 x 10~ 15 a tr cosh 4/3 [27.09 + 2.315 x 10~ 43 (t - t tr )\ (656) 

using (50a) and (64). 

(65b) shows acceleration for t > t tr . It is interesting to see a(t), given 
by (65b), will reduce to de-Sitter like expansion asymptotically. 

Further using (50a) and a = a(t ) in (65b), it is obtained that 

t tr ~0.183t - (65c) 

4. Summary Here, cosmology is obtained from the variable gravita- 
tional constant and its fluctuations around the stable state G = Gn- 
Variation in G is manifested by a scalar <p{x) given as 
The theory is based on <p and the Ricci scalar curvature R only. So, 
the action contains Einstein-Hilbert like term with G depending on 
as well as kinetic energy and self-interaction potential term for it. The 
model of the universe begins at the Planck scale being the fundamental 
scale. 

This model provides an interesting cosmology. It is found that trace 
of the gravitational field equations yields the Friedmann equation with 
dark matter and a fluid term behaving as quintessence dark energy. 
These terms emerge as combined effect of and R. Here, it is in- 
vestigated that the very early universe inflates for an extremely small 
period of time being 5.7 x 10~ 18 sec and exits from the inflationary 
phase at t c — tp + 5.7 x 10~ 18 sec(t p = Mp 1 is the Planck time), when 
4> = 4> c = a/3/87t x 10~ 12 Mp. A phase transition takes place at this 
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epoch and a huge amount of energy is released as latent heat ther- 
malizing the universe upto the temperature 2.45 x 10 14 GeV After the 
inflationary phase, universe decelerates driven by emitted radiation. 
But, due to the effect of variable G, a(t) ~ t 1 / 3 , which is not consistent 
with standard radiation-model. 

After sufficient expansion, when temperature of the emitted radi- 
ation red-shifts to the temperature 246GeV (the electro- weak scale), 
another phase transition takes place when rolls down to the stable 
state = a/3/87t|Mp|. Major portion of the released energy reheats 
the universe slightly and the temperature is raised by 5.22 x 10 5 x 
and an extremely small portion of the released energy of the order of 
2.95 x 10~ 52 GeV 4 is converted to dark energy caused by fluctuations 
around the stable state = y/3/8n\Mp\. Thus, after the second phase 
transition, G is no longer a variable and acquires the constant value 
G = Gn- Radiation still dominates the cosmic dynamics. But, in 
the post-second phase transition era, radiation-dominated phase de- 
celerates ~ t 1 ! 2 being consistent with the standard model. Thus, the 
standard model of cosmology is recovered at the epoch of the second 
phase transition and the radiation emitted at this epoch is identi- 
fied as the cosmic background radiation having the present temper- 
ature 2.73°^. Thus, in this model, standard cosmology begins with 
temperature 246GeV + 5.22 x lO 50 -^. It happens so at the red-shift 
z = 4.47 x 10 15 . 

As mentioned above, fluctuates around the stable state = ^/3/8tt 
x|Mp| after the second phase transition, which oscillates for a very 
small period till a ~ 1.004a*. It is found that phantom dark energy 
is seeded by the released energy during fluctuation, which grows later 
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(after fluctuation ends) in the expanding universe such that the present 
dark energy density is 171.137 times the phantom dark energy at a ~ 
1.004a*. This phenomenon takes place during radiation-dominance. 

As radiation density falls faster than the dark matter density dark 
matter dominates the radiation density at red-shift ~ 1091.98. It is 
interesting to see that this red-shift is very closed to the red-shift 109812 
for decoupling of matter from radiation as per WMAP observations 
[25]. During the dark matter-dominance, the universe decelerates as 
~ t 2//3 , which is consistent with standard matter-dominated model. 

It is discussed above that the phantom dark energy grows during 
radiation-dominance and matter-dominance. The dark matter density 
decreases in the expanding universe. So, obviously on sufficient growth 
of a(t), dominance of phantom over dark matter is expected. In the 
present scenario, it is found that phantom overpowers dark matter at 
red-shift z tr = 0.46, which is within the range 0.46 ± 0.13 obtained by 
16 Type Supernovae observations [2]. This transition takes place at 
time ttr = 0.13to- At this epoch, phantom energy with growing density 
in expanding universe gives a jerk. As a consequence, the decelerating 
universe is accelerated. 

Interestingly, EOS parameter w p hde for phantom energy, created by 
fluctuation of 0, depends on the scale factor a(t). It is obtained that 
Wphde — > —1 and Pphde — > l-9747pp Me as a — > oo. Also, a(t) reduces to 
the de-Sitter form asymptotically. These results show that ultimately 
phantom energy will evaporate contributing to the vacuum energy. 

The following are salient features of this model, 
(i) Quintessence dark energy in the early universe, visible radiation, 
dark matter and phantom dark energy are caused by the scalar field 
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4> = ^/3/87rG _1 / 2 . Moreover, phantom dark energy is created due to 
fluctuations of around the stable state = ^/3/87r|Mp|. 

(ii) Two phase transitions take place in this cosmological picture. The 
first one takes place at the GUT scale when the early universe comes 
out of the inflationary phase and the second one takes place at electro- 
weak scale. This epoch heralds the standard model of the universe with 
background radiation having temperature ~ 246GeV with its remnant 
having temperature 2.73°-K" in the present universe. 

(iii) Here, it is shown that phantom dark energy is created by fluctu- 
ations in the stable state = y/3/8nMp = . Thus, this 
model provides a possible way of creation of phantom energy almost 
in the beginning of the mtter-dominated phase of the universe. It is 
found that growth of phantom density is very fast during radiation- 
dominance in the post-second phase transition era, but it slows down 
later on. It is found that phantom density will become almost double 
of the present dark energy density, when a — > oo. It is interesting 
to obtain that phantom density will evaporate and will contribute to 
vacuum energy ultimately. 

(iv) The scale factor a(t) is free from big-rip problem. 

Thus, this model yields a physically viable and unified cosmological 
picture from early universe to the future universe. Results are consistent 
with observations. Most interestingly, a possible way for creation of 
phantom dark energy is suggested here. 
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